AN INTERESTING CLASS OF 
PARTIAL DIFFERENTIAL EQUATIONS 
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WEN- AN YONG 

Abstract. This paper presents an observation that under reasonable conditions, many 
partial differential equations from mathematical physics possess three structural proper- 
ties. One of them can be understand as a variant of the celebrated Onsager reciprocal 
relation in Modern Thermodynamics. It displays a direct relation of irreversible processes 
to the entropy change. We show that the properties imply various entropy dissipation 
conditions for hyperbolic relaxation problems. As an application of the observation, we 
propose an approximation method to solve relaxation problems. Moreover, the obser- 
vation is interpreted physically and verified with eight (sets of) systems from different 
fields. 
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1. Introduction 

The goal of this paper is to draw attention to a class of partial differential equations 
(PDEs) of the form 



;i-i) 



R"^ X 



Here U is the unknown n- vector- valued function of {x,t) = {xi,X2,--- ,Xd,t) G 
[0, -|-oo), taking values in an open subset G of R" (called state space) ; Q{U) and Fj{U){j = 
1,2, ■■ ■ ,d) are given n-vector-valued smooth functions of U & G; and the subscripts t 
and Xj refer to the partial derivatives with respect to t and Xj, respectively. 

As fundamental PDEs and as intermediate models [H |T51 [12] between the Boltzmann 
equation [2J and hyperbolic conservation laws [5], systems of first-order PDEs with source 
terms describe various irreversible processes of scalar type |llj . Important examples occur 
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in chemically reactive flows [9] , radiation hydrodynamics [T71 [22] , invisicid gas dynamics 
with relaxation [31], nonlinear optics [12], and so on. 

Since the last decade, PDEs of the form (11. ip have attracted much attention. See [20 tl28| 
[25] and references cited therein. One of the main interests is to identify a set of structural 
properties (axioms) that are satisfied by most of important equations from applications 
and, meanwhile, provide a convenient framework for the development of mathematical 
theories. In this regard, two stability conditions and various entropy dissipation conditions 
have been proposed in [26] and [21 [191 [29l [131 [25] , respectively. See also [281 [23]. All those 
conditions are generalizations of the well-known subcharacteristic condition [T6] for (11.11) 
with n = 2 and d = 1. For (II. ip . such a condition is the same in spirit as the H-theorem 
for the Boltzmann equation [2] and as the entropy condition for conservation laws [1]. 

In this paper, we present an observation that under reasonable assumptions, many 
equations of the form (II. ID from mathematical physics fall within a class characterized 
with the following three properties. (I) Every system in the class admits a strictly convex 
entropy function [TUl [S] , (H) the source term can be written as a product of a non-positive 
symmetric matrix and the corresponding entropy variable, and (III) the symmetric matrix 
has a constant null-space. 

The first property is the well-known entropy condition for conservation laws and cor- 
responds to the classical principles of thermodynamics. Property (II) can be understand 
as a variant of the celebrated Onsager reciprocal relation in Modern Thermodynamics 
[TTl [Hj and implies the second law of thermodynamics. It displays a direct relation of 
irreversible processes to the entropy change. Property (III) expresses the fact that physi- 
cal laws of conservation hold true, no matter what state the underlying thermodynamical 
system is in (equilibrium, non-equilibrium, and so on). 

We will verify the three properties for eight (sets of) systems of the form (II. ip arising 
in gas dynamics with damping or with relaxation, nonlinear optics, radiation hydrody- 
namics, chemical reactions, kinetic theories (both moment closure systems and discrete 
velocity models), and so on. Furthermore, we show that the properties ensure a uniquely 
defined Maxwellian and imply various entropy dissipation conditions in the literature for 
hyperbolic relaxation problems. Thus, all the general results in [27 1[30l[T8l[23l[29l[T3l [25l[5] 
apply to the aforementioned fields. 

We notice that the examples in Sections 8 and 9 have a common kinetic origin. More- 
over, we know from [S] that the chemical systems in Section 7 have a similar origin. It 
would be interesting to include the radiative gas example of Section 6 in the same basket, 
with the idea that radiation is associated with particles (photons) collisions. In other 
words, we show that the kinetic theory yields the Onsager relation for gas mixtures. 

As an application of our observation, we propose an approximation method to solve the 
relaxation problems. The accuracy of the method is analysed for initial value problems 
with smooth initial data, by using the results in [27]. In this analysis, an important 
ingredient is a continuation principle for hyperbolic singular limit problems (Lemma 9.1 
in ^28j and the appendix in [Ij). Further analysis and applications of the approximation 
method are desirable. 

Our discussions indicate that the above three properties have a solid basis, from both 
mathematical and physical points of view. Thus, it seems reasonable to take the properties 
as requirements in construction of new mathematical models for irreversible phenomena. 
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The paper is organized as follows. In Section 2 we present the observation and discuss 
its mathematical consequences. Section 3 is devoted to the approximation method. Some 
physical interpretations are given in Section 4. The rest of the paper deals with the eight 
(sets of) examples. 

2. An Observation 

The main contribution of this paper is the following observation that under reason- 
able assumptions, many partial differential equations (PDEs) of the form (11.11) arising in 
mathematical physics admit the following structure: 

(I) . There is a strictly convex smooth function r]{U) such that riuu{U)Fjjj{U) is sym- 
metric for all f/ G G and all j. 

(II) . There is a symmetric and non- negative definite matrix C = C{U) such that 

QiU) = -C{U)MU). 

(III) . The null space of C{U) is independent of U EG. 

Throughout this paper, riu{U) should be understand as a column vector. We will use * 
as a superscript to denote the transpose operator. 

Recall that (I) is the classical observation due to Godunov and Friedrichs and Lax 
[6] for conservation laws 

d 

i=i 

In what follows, we point out several important conclusions of the above observation. 
First of all, we show 

Proposition 2.1. The observation is invariant under linear transformations of the form 

V = PU. 

Here P is a constant and invertible n x n-matrix. 

Proof. Let rj = ri{U) be the strictly convex function in the observation. We need to show 
that the observation holds true with fj{V) = ri{P^^V) for 

(2.1) Vt + ^(PF,(P-V)).^. = PQ(P"V). 

j 

To this end, we compute 

7]u = P*fjv and r/uu = P*VvvP- 

The second equality indicates that fjvviV) is positive definite, since so is rjjjjjiU). The 
latter is equivalent to the strict convexity of f]{U). Therefore, fj{V) is strictly convex. 
Since 

fivv{PF,{P-'V))v = {P-*vuuP-'){PF,uP-') = P-*VuuFjuP'' 
is symmetric, (I) follows. Moreover, since 

PQ{P-W) = -PCr]u = -PCP*f]v 
and P is independent of U, (II) and (III) follow. This completes the proof. □ 
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Let r be such a constant that (n — r) is the dimension of the null space in (III). In 
Proposition 12. ![ we take P to be such a matrix that its first {n — r) rows span the null 
space. Then the next proposition becomes obvious. 

Proposition 2.2. Under the observation, there is a constant invertihle matrix P and a 
symmetric positive definite r x r -matrix 1{U) such that 

PC{U)P* = (im^(0(„_,)x(„-o, ^(U)). 

Here and below we denote by Ox the origin of H.^ to avoid possible confusions. 

With P from Proposition 12.21 the source term in (12.11) obviously has the form 

FQ(F-'V) = (»-) 

with q{V) G R''. Accordingly, we introduce the partition 

V = 





and rewrite (12. ip as 
(2.2) 

Notice that 

(2.3) q{u, v) = -l{u, v)ri^{u, v). 

Because P is invertible and constant, (12.21) with (12.31) is equivalent to (II. ip . Thus, 
one may assume that (11.11) is already in the form (12.20 with (12.31) . In particular, the 
symmetric matrix C in (II) has the block-diagonal form in Proposition 12.21 and 1{U) 
therein is symmetric and positive definite. 

With the equivalent form (12.21) . we turn to the following theorem, which is related to 
the so-called equilibrium manifold 

S := {UeG:Q(U) =0}. 

Theorem 2.3. Assume that U.l\) possess the three observed properties and the state space 
G is convex. Then for every U (z G, there is at most one point in £, say M = M{U), 
such that U — M{U) is in the orthogonal complement of the null space. Moreover, there 
are two positive functions c{U) < C{U), defined in the domain of M = M{U), such that 

(2.4) c{U)\U - M{U)\ < \Q{U)\ < G{U)\U - M{U)\. 

Here \X\ denotes the Euclid norm of matrix or vector X . 

Note that M{U) may be defined only in part of the state space. A prototype of MU) 
is the Maxwellian in kinetic theory [2]. 

Proof. Assume that for a certain U & G, there are two points Mi, M2 G S such that both 
U — Ml and U — M2 are in the orthogonal complement of the null space. Let P be the 
invertible matrix in Proposition 12.21 and set 

PU = {u, v), PMi = (Ui, Vi), PM2 = iu2, V2). 
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Then we have q{ui,Vi) = = q{u2,V2) and ui = u = U2. Since l{u,v) is positive 
definite, from (12.31) it follows that q{u,v) = if and only if fj^{u,v) = 0. Thus we have 
fjj,{u,vi) = = fjy{u,V2)- This contradicts the strict convexity of fj and thereby proves 
the first part of the theorem. 

For the inequalities in fl2.4p . we denote by m{U) the last r components of M{U). Since 
fj = fj{U) is strictly convex, there is a symmetric and positive definite r x r-matrix a{U) 
such that 

UPU) = UPU) - UPM{U)) = a{U){v - m{U)). 

Thus, from ([23]) we have q{PU) = -l{PU)a{U){v -m{U)). Now the inequalities in fl^ 
can be easily verified with 

c{U) ={\P\\p-'\\l-\PU)a-\U)\)-\ 
C{U) =\P\\p-^\\l{PU)a{U)\. 
This completes the proof. □ 

The following theorem provides connections of the observation to the existing entropy 
dissipation conditions in the literature for hyperbolic relaxation problems. 

Theorem 2.4. Assume that ( fi.ij) possess the three observed properties. Then the follow- 
ing conclusions hold: 

(1) 7]Ij{U)Q{U) < -X-\U)\Q{U)\^ for all U e G, where X{U) is the maximum 
eigenvaue of C{U). 

(2) Q{U) = if and only if ri^{U)Q{U) = if and only if f]u{U) is in the null space. 

(3) 7]lj{Ue)Q{U) = for any U^UeEG with satisfying Q{Ue) = 0. 

(4) For U satisfying Q{U) = 0, Qu{U)ri^]j{U) is symmetric and non-positive definite, 
and its null space coincides with that of C{U). 

(5) For the equivalent version Ii2.2^) with Ii2.3\) . r x r-matrix q^{u,v) is invertible for 
{u, v) satisfying q{u, v) = 0. 

Proof. (1). Since C{U) is symmetric and non-negative definite, it is clear that C^{U) < 
\{U)C{U). Therefore we have 

V*u{U)Q{U) = -vh{U)C{U)vu{U) < -X-'\C{U)vum' = -X-\U)\Q{U)\\ 

Note that the symmetry of C{U) is crucial to the inequality. 

(2) . From (1) it follows immediately that Q{U) = if and only if r]lj{U)Q{U) = 0. Since 
t]Ij{U)Q{U) = -T]lj{U)C{U)r]u{U) and C{U) is symmetric, t]*jj{U)Q{U) = is obviously 
equivalent to that rjuiU) is in the null space of C{U). 

(3) . It follows from (II) that for f/g satisfying Q{Ue) = 0, rjuiUe) is in the null space of 
C{Ue). Thanks to (HI), TjjjiUe) is in the null space of C{U) for any U E G. Thus, we have 

(2.5) C{U)vu{Ue) = 

and thereby T]lj{Ue)Q{U) = for any U,Ue e G with Ue satisfying Q{Ue) =0. 
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(4) . From (II) and ([23]) it follows that 

(2.6) Q{U) = -C{U){vu{U) - VuiUe)). 

With this relation, it is clear that 

Thus (4) becomes obvious. 

(5) . Recall the block- diagonal form of C{u,v) for (12. 2p with (12. 3p . It follows from (2) 
that fiy{u,v) = for {u,v) satisfying q{u,v) = 0. Thus, we compute from (12. Sp that at 
{u,v) satisfying q{u,v) = 0, qy{u,v) = —l{u,v)fjyv{u,v). Note that fjyy{u,v) is positive 
definite, since ri{u, v) is strictly convex. This, together with the positive definiteness of 
l{u, v), implies the invertibility of g„(u, v). Hence the proof is complete. □ 

By Theorem 12.41 if a system of PDEs possesses the three observed properties, then it 
satisfies all the entropy dissipation conditions in [3l |19l |29l |13l |25] . In fact, all the existing 
conditions consist of (I) and some additional requirements. They are rj^{U)Q{U) < 0, 
which is implied by (1) of Theorem 12. 4[ and (2) of Theorem 12.41 in [3]. In [19], there is only 
one additional requirement which is ri^{U)Q{U) < — the second law of thermodynamics. 
In [21], the additional requirements are (1) and (5) of Theorem 12. 4^ while they are (2) 
and (4) in [T3|. The entropy dissipation condition in [23] is that in [22] with |(5(f/)| in (1) 
replaced by \U — M{U)\. For this see the inequalities in (12. 4p . Moreover, it was shown 
in [29] that the entropy dissipation condition therein implies the stability conditions in 
[261 [27]. Thus, all the general results in [3 [131 [HI [231 [25l [271 [291 [30] are vahd for PDEs 
of the form (II. ID which possess the three observed properties. 



3. An Approximation Method 
Consider (II. ip with a small parameter e > 0: 

(3.1) Ut + Y,Fm.,=QiU)/e. 

j 

This is the so-called relaxation problem. Assume the three observed properties in the 
previous section hold for this relaxation system (13. ip . Then the scaled system can be 
rewritten as 

j 

As an application of our observation, we show in this section that the above relaxation 
system can be approximated with the following simplified system 

(3.2) Ut + X]F,(t/)., = -£iU,)vuiU)/e 

3 

as e tends to zero. Here f/* G G is arbitrarily fixed. 
To this end, we use the equivalent form (12. 2p : 



(3.3) 1 + ^ 
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Here 

q{u,v) = -l{u,v)ri^{u,v) 

as in (12.31) . It is not difficult to see that as e goes to zero, the formal limit of solutions to 
(13. 3p solves the following so-called equilibrium system 

ut + ^fj{u,v)^.=0, q{u,v) = 0. 
j 

This system consists of differential and algebraic equations. 

Since q{u,v) = —l{u,v)fj^{u,v) and l{u,v) is positive definite, the equilibrium system 
is equivalent to 

j 

This system is independent of l{u,v). By Theorem 12.31 the algebraic equations define v 
as a unique function of u, say, v = h{u). Here we assume that the domain of h{u) is 
non-empty and open! Thus, the equilibrium system becomes 

(3.4) + ^/j(m,/i(u))^^ = 0, v = h{u). 

j 

It is remarkable that h{u), and thereby the equilibrium system, is independent of l{u,v)\ 
As is pointed out in the previous section, relaxation system (13.11) satisfies the stability 
conditions in [261 [27], for it possesses the three observed properties. Thus, Theorems 6.1 
and 6.2 in [27] apply here: For smooth initial data, there is a finite and e-independent time 
interval [0, T] such that the initial value problem of (13.11) has a unique smooth solution 
U" = U''{x,t) defined for t G [0,T] and satisfying 

(3.5) ^--^"(mJ^^W 

in a certain Sobolev space, as e goes to zero. See [21] for details. Here u solves the 
corresponding initial value problem of the equilibrium system in (13. 4p . In addition, we 
have assumed for simplicity that the initial data take values in equilibrium and thereby 
initial- layers do not appear. Recall from [27] that the time interval [0,T] is the life-span 
of the smooth solution u. 

Note that (13.41) is also the equilibrium system for the corresponding equivalent version 
(12. 2p of the simplified system (13. 2p . The latter obviously possesses the three observed 
properties. Thus, we see that with the same initial data, the simplified system has a unique 
smooth solution Lf^ defined in the same time interval and having the same expansion 



(3.6) U' = p-U +0(e 




as e goes to zero. Here the key points are the same equilibrium system and the same time 
interval. The latter is attributed to a continuation principle for hyperbolic singular limit 
problems (Lemma 9.1 in [28j, see also the appendix in [T]). Consequently, we see from 
and dSSD that 

U'-U' = 0(e), 
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in a certain Sobolev space, as e goes to zero. 

In conclusion, we have shown that for small e, relaxation systems (13. ip and (13.21) are 
close to each other in a finite and e-independent time interval for initial value problems 
with smooth data. The above discussion suggests an approximation method to solve the 
original relaxation system (13. ip . Further analysis and applications of this approximation 
method are desirable. In particular, it would be interesting to study the closeness for 
specific systems in the regime of non-smooth solutions. 

4. Physical Interpretations 

In this section, we give some physical interpretations of the three observed properties in 
Section 2. Recall that for a thermodynamic system inside which n irreversible processes 
occur, the infinitesimal entropy change ^5" due to the processes can be expressed as a sum 
of two parts: 

dS = deS + diS. 

Here dgS is the part supplied to the system by its surroundings, and diS is that produced 
inside the system. It is well known (see, e.g., [TT]) that deS corresponds to the fiux terms 
in (11.11) and diS to the source term. The second law of thermodynamics states that diS 
is zero for reversible processes and positive for irreversible ones. 
Based on our observation, (II. ip can be rewritten as 

d 

(4.1) U, + J2f,{U),^ = -C{U)vu{U). 

This form relates irreversible processes directly to the entropy change rju- 

Recall that the physical entropy S is equal to —r] and its existence is guaranteed by the 
classical principles of thermodynamics pTl[H] . This explains why the classical observation 
(I) has a solid basis in thermodynamics. The Gibbs relation on the total differential of t] 
(or S), in this general level, reads as 

(4.2) d7] = r]u{U) ■ dU, 

where the dot " ■ " between two vectors means the scalar product. The usual Gibbs 
relation 

edS = de + pd(-) H 

P 

is a slight rearrangement of (14. 2p . Here 6 is the temperature, e is the specific internal 
energy, p is the pressure, p is the density, and the dots come from other possible internal 
variables. 

Property (II) very much looks like the celebrated Onsager reciprocal relations in Non- 
equilibrium Thermodynamics [HI [T3] , if one understands the source terms as irreversible 
fiuxes and the entropy variables as thermodynamic forces or affinities. However, it is 
slightly different from the Onsager relation. Firstly, it seems new to choose the entropy 
variables, instead of their linear combinations, as thermodynamic forces. Secondly, unlike 
the Onsager relation, 

QiU) = -C{U)r,u{U) 

is a nonlinear relation between Q{U) and rjuiU). In fact, the matrix C = C{U) depends on 
U. Because the entropy function is strictly convex, there is a one-to-one correspondence 
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between U and the entropy variable rjuiU) (see |E] for a proof of this fact). Thus, C 
depends on the entropy variable, which plays the role of affinities here. 

By the way, it is well known (see, e.g., [21], page 125-126) that there are difficulties in 
choosing the thermodynamic forces and fluxes when applying the notion. Here we have 
proposed an unconventional but unambiguous choice of the couple. 

Furthermore, we recall (12. 6p and deduce that for any U,Ue G G with Uf. satisfying 

g([/e) = 0, 

Q{U)=-C{U){Tiu{U)~Tlu{U:)) 

= - C(Ue){r,u{U) - 7]u{Ue)) " {C{U) - C{Ue)){Vu{U) - VuiUe)) 

= -C{Ue)vu{U) + 0{\U-Ue\^). 
Neglecting the higher-order term, we obtain a linear relation 

Q{U) = -CiUMU) 

between Q{U) and riu{U). Because C{Ue) is symmetric, this is the Onsager reciprocal 
relation if one considers the source terms as irreversible fluxes and the entropy variables 
as affinities. 

As to Property (III), we recall the equivalent form 02.21) of fll.lj) . In fl2.2p . the first (n—r) 
equations represent (n — r) conservation laws. Note that r might not have been a constant 
without assuming (III). In other words. Property (III) expresses the fact that the physical 
laws of conservation hold true, no matter what state the underlying thermodynamical 
system is in (equilibrium, non-equilibrium, and so on). 

5. Four Specific Examples 

From this section on, we will verify the three observed properties in Section 2 for a 
number of systems of the form (11. ip arising in applications. This section contains four 
comparatively simple examples. 

Example 1. Mult i- dimensional Euler equations of gas dynamics with damping: 

Pt + div(pM) =0, 
{pu)t + div(pM m) + V]9(p) = — pu. 
As usual, p = p{x, t) stands for the density and u = u{x, t) is the velocity. This system is 
of the form with U = (^^^^ G R'^+^ 
It is well known that function 



■n{U) = ^-^+ J J ——dadT. 

is a strictly convex entropy for the above system in the classical sense (I). By computing 
riu{U), we see that 

Q{U) = -dia.g{0,pI,)7]u{U), 

where Ik is the unit matrix of order k. Thus, the properties (II) and (III) obviously hold 
with C{U) = diag(0,p/d) for p > 0. 
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Next three examples all have the form (12. 2p with r = 1. For such a system, if there is 
a function r} = r]{u, v) satisfying Property (I), then the observation is obviously true with 

jC{U) = - ' ^ diag(0(n_r) x(n-r), 1) 

{Ox is the origin of R^), provided that 

5.1 — -<0 

Vv{u,v) 

for all {u, v) under consideration. The inequality (15.11) is a stability condition for the 
corresponding systems. 

Example 2. A 3-D quasilinear system for nonlinear optics: 

A - V X = 0, 
Bt + VxE = 0, 

Xt = \E\''-x 

with D = {1 + x)E- See [12J for an explanation of the equations above. The state space 

here is G = {{D,B,x) : I) G R^ 5 e R^, x > O} C R^ 
Set 

U= (b 

In [12], Hanouzet and Huynh showed that function 

r^{U) = {l + xr'\D\' + \B\' + xV2 

is a strictly convex entropy in the classical sense (I) in order to study the corresponding 
relaxation limit of the above system. By computing riu{U), we see that 

Q(f/) = -diag(06x6,l)w(^^)• 
Thus, the observation is true with C{U) = diag(06x6! !)• 

Example 3. 1-D Euler equations of gas dynamics in vibrational non-equilibrium (in 
Lagrangian coordinates) : 

ut - =0, 
Ut + Px =0, 

(e + y)* + {pu)x =0, 

qt =uj{9i) - uj{92). 

See [31] for an explanation of the equations above. 
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For this system, we know from [31] that there is a strictly convex function rj = //(z/, u, e+ 
q) such that Property (I) holds and 



Then we have 



uj(ei) -00(62) 
QiU) = -9,92 ^ J / 'M iag(03x3,l)w(^)- 

(71 — (72 



Thus, the observation is true with 

m) = 9,92 ^ ' / M iag(03x3, 1), 

(71 — (72 

for uj = uj{9) is strictly increasing |31j . 

Example 4. 1-D Euler equations for isothermal motions of a viscoelastic material (in 
Lagrangian coordinates): 

ut - =0, 
Ut + Px =0, 
{p + Ep)t = -p- g{p). 

See [24] for an explanation of the equations above. 
For this system, we know from [2^ that function 

/—p—Ev 
h-\a)da, 

is a strictly convex entropy in the classical sense (I). Here is the inverse of /i(z/) = 
g^u) — Eu, which exists under the so-called sub characteristic condition 

(5.2) < g^iu) < E. 

Since 

r/p(t/) = h-\-p-Ev)-p, 

we have 

= -FT(^J^Zli^diag(0...,l).„(C'). 

Thus, the observation is true with 

= /.-H-i>-''it)-i-HM^)) ^^'^^°^^^''^' 

for h{u) = g{u) —Eu is strictly decreasing [21] under the subcharacteristic condition (15.21) . 
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6. Radiation Hydrodynamics 

In this section, we consider discrete-ordinate models of the Euler equations for radiation 
hydrodynamics [I71[22], which are of the form (11.11) with 



U 



pVi 
PV2 

pE 

h 



F,{U) 



( \ 

pviVj + 5ijp 

pV2Vj + 62jP 

pvsVj + Ssjp 
pEvj + pvj 

p]h 



\ 



Q{U) 



J 



( \ 






CpYi=^{ii-B{Qy) 

-Pih - 3(9)) 
V -pih-BiO)) J 



Here p is the density, vj is the velocity in the j^^ direction, E = e + |f P/2 with e the 
specific internal energy, Ii is the radiation intensity in the direction p'- = (/i'^, /X2, /^a), 
p = p{p, e) is the pressure, 6ij is the standard Kronecker delta, C is a positive constant, 
and B = B{6) is the Planck function of temperature 6. 

For this system, the state space is (0, oo) x R'^ x (0, 00)^+-'^. Since the basic assumptions 
of radiation hydrodynamics are not valid at low temperatures, we restrict the temperature 
domain to [6q, 00) with 6*0 > a constant. 

Recall that B = B{0) > is strictly increasing with respect to > 6*0. We denote by 
h = h{y) the inverse function of B{6), that is. 



(6.1) 



e = h{B{e)), 



y e> On. 



Note that b = b{y) is strictly increasing. Moreover, it is smooth if so is -B = B{6). 
Define 



(6.2) 



^{U) = -ps{p,e)-Cj2 



1=1 



B{eo) 



dy 
Ky) 



with s = s{p,e) the specific entropy. It is straightforward to verify that this rj is strictly 
convex. Since the system is the classical Euler equations coupled weakly to L linear 
transport equations, t] is obviously an entropy function for the system. Namely, Property 
(I) is verified. 
Note that 



Vpe{U) 



1 

1' 



and set 



B{9) 



C 



e-^-h-^ii)- 



Then it is not difficult to see that 



QiU) = -C{U)r,u{U) 
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with 



(6.3) 



jC{U) = p 



X 4 


1 


1 


O4 X 1 


04X 1 


O4X1 \ 


Olx4 






-02 


-CTg ■ 


-ctl 


Olx4 


-CTi 













Olx4 


-0-2 













Olx4 


-CT3 













V0lx4 


-0-L 














Since > for all Z, this C{U) is symmetric and non-negative. Moreover, its null space 
is 



span 



61,62,63,64,65 +cy^ei 



/>6 



Conse- 



which is independent of U. Here 6^ is the k column of the unit matrix I(l+5)- 
quently, the properties (II) and (III) are also verified. 

7. Chemically Reactive Flows 

Most of this section is taken from except the verification of the properties II and 
(III). For multi-component reactive flows, if we neglect external forces, diffusion of mass, 
heat conduction and viscosity, and but retain the chemical reactions, the flows are de- 
scribed with PDEs of form (11.11) . where 

/ Pi \ / PiVj \ ( miui \ 



(7.1) 



U 



P2 

pvi 

PV2 
PV3 

\pE/ 



PusVj 

pviVj + 6ijp 

PV2Vj + S2jP 

pvsVj + hjP 
V pEvj + pvj ) 



Q{U) 



ni2UJ2 








Here pk is the density of the k species, Ug is the number of the species, p = pk is 
the total density, Vj is the mass averaged flow velocity in the j^^ direction, E = e + \v\^ /2 
with 6 the specific internal energy of the mixture, p is the pressure, 5ij is the standard 
Kronecker delta, is the molar mass of the k^^ species (known constants), and ojk is the 
molar production rate of the k^^ species. The system of will be closed by specifying 6,p 
and ujk as functions of the natural variable 

(7.2) Y = (pi,P2, ■ ■ ■ ,pns,Vl,V2,V3,ey 

with 9 the absolute temperature. 

We will specify Uk later, p and e are given as in [9] . For p, we denote by Rg the universal 
gas constant, write = Rg/ruk and then define 

(7.3) p = Oy^jkPk- 
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e is taken as the weighted average of the specific internal energy of the k^^ species: 

(7.4) pe=^Pfcefc, 

k 

where 

(7.5) efc =e° + / c^k{y)dy. 

Here is the specific internal energy of the k^^ species at the reference temperature 
^0 > 0, and c^k = c^kiO) are given smooth functions of 6^ G [6^0, oo), denoting the specific 
heat at constant volume of the k^'^ species and satisfying YmYikfi{cvk{0)} > 0. 

The state space for the natural variable Y defined in (17. 2p is (0, cxd)"" x R'^ x [9o, oo). 
For the conserved variable U defined in (17. ip . it is 

G = {Ue R"^+^ : [/fc > for 1 < A; < and f/„,+4 > 0(t/i, f/2, ■ ■ ■ , t/n.+s)} , 

where 



Since is a convex function, the state space G is convex. 
Introduce 

(7.6) Skipk,0) = sl+ [ ^"^^^^ rfy - rfc In (— 



do y 



where s° is a constant, and define 

(7.7) v{U) = -^pkSk{pk,i 



k 

We show that this 77 = ri{U) is an entropy function in the classical sense (I). Since 

12 



pE = ^Pk (^4 + ^ c^k{y)dy^ + 



p\v\ 



due to (17.41) and (17. 5p . we compute to obtain 



Ou = ^^^PkCvkiO)^ 



2 ei, — - e„^, -t;2, -^^3, 1 ) • 



Thus, it follows from (17. 7p and (17.60 that 

(7.8) VU = -^{Pl- — , ■ ■ ■ , /^n. - f 1, W3, -1 

where 

(7.9) Pk = ek + TkO - SkO 
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denotes the chemical potential of the k^^ species. On the other hand, by the definitions 
of U and Y in (17. ip and (17.21) . we compute 



dY 

dU 



Ins Oris X 3 OjisXl 

-p-it;(l,l,--- ,1) p-1/3 O3XI 
0u 



with V = {vi,V2,Vs)* , and thereby 
drju drju dY 



rjuuiU) 



dU dY dU 

Osxns 
Olxns 



0lx3 
Oris X 3 



OrisXl 

Oaxi 



0lx3 J2k PkCvkiO) ^ 



dY 
dU 



dY 
dU 

> 0. 



Therefore, ri{U) is strictly convex. Moreover, since 

/OnsXl 



due to (17.11) . we have 



P2 

Pus 

pv 

\pE + pj 



du{v ■ 




Therefore, it follows from ([LHD, (Q, (II3D, dUD and ([TJD that 

Hence, rjlU) is a strictly convex entropy function for the system in the classical sense (I). 

Next we turn to specifying uj^ by following [S] again. Let the system have reversible 
reactions for species: 



for i = 1,2, ■■■ ,72^. Here Sk is the chemical symbol for the /c*^ species, and and 
i/^'j are the stoichiometric coefficients of the /c*^ species in the i^^ reaction. The molar 
production rates Uk are the Maxwellian production rates obtained in the kinetic framework 
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of the "slow reaction regime" or in the "tempered reaction regime", when the chemical 
charactersitic times are larger than the mean free times of molecules: 



(7.10) U = {Ui, ■ ■ ■ , iOnsT = X] Tii'^li, ^2i, " " " , J^risi)* = 



Here = i^'L — v',- and is the rate of progress of the i*'* reaction: 



fc ^ '^^ k ^ " 



ki 



where KfiiO) and KniO) are the direct and reverse constants of the i^^ reaction, respec- 
tively; and 



K,0) :=exp i-Y^{rue)-h 

\ k 



(7.12) = K,0) := exp I - > ;(r,^^)-Vfci/Xfc(m,, 9) 



with fik{rnk,0) the chemical potential (17. 9p at the unit concentration: pk/riik = 1. 
It is well know that 

^mfcCUfc = 0. 

k 

In fact, let ne be the number of elements involved in the system and denote by en the 
number of the element in the fc*^ species. We have the element conservation relations 

(7-13) ^^'ki^ki = ^i^ki^ki 

k k 

for i = 1,2, ■ ■ ■ ,nr and / = 1, 2, ■ ■ ■ , ng. On the other hand, the species molar mass is 
related to the elemental masses ai by the relation 



I 

Hence it follows from (17.10p and (I7.13P that 

^ rrikUJk = X] aiekiTiUki = ^ am ^ enii^ki - v'ki) = 0. 

Similarly, we have 

(7.14) ^rnfci/fc^ = 0. 

k 

To see the properties (II) and (III), we set 

M = diag(mi,m2, ■ ■ ■ ,mnj, y = {Rg9)~^{ni, H2, ■ ■ ■ , fJ^nJ' 

and 

A, = Kf,{9) n (-^) [ exp(a < 3^, Mz/, >)da > 0. 
, \fnk/ Jo 
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It follows from fmTD . fl7l2|) . (\7^ . ([71]), fITJOD and (04]) that 
r, = - Aiiy:My 

/' D A ,,* A/f 



Note that the reaction rates Tj depend exponentially on the chemical affinities i>*My, due 
to 

exp{<y,Mu,>) - I 



/ exp(o- < 3^, Mui >)da 
Jo 



Moreover, we set 

V = {vi,V2,--- and A = diag(Ai, A2, ■ • ■ , A„J. 

Then we deduce from (O), flTrTOD and (EHD that 

Q{U) = -R;'dmg{MVAV*M,0,^,)vuiU). 

Since MV is a constant matrix and A is positive definite, the null space of MV AV* M is 
independent of U. Hence the observation is verified with 

C{U) = R-^diag{MVAV*M, 04x4)- 

Finally, let us mention that some statements of Proposition 2.4 are also discussed in 

8. Moment Closure Systems 

Moment closure systems in kinetic theories are PDEs of the form (11 .11) . In this section 
we show that our observation holds for the exponentially based closure systems in [T5] 
corresponding to the Boltzmann equation 

(8.1) /, + ^.V./= / {fj:-ff')B{u;,^,Odujd^'. 

Here / = f{x,t,C,) > denotes the kinetic density of particles at the position-time- 
velocity point (x, t,C,) E IV^ X R_|_ x R*^, the dot " ■ " between two vectors means the scalar 
product, = /(x,t,^^),/^ = fix,t,Q and /' = /(x,t,^') with 

e. = e-^-(e-e'v and e: = r+^-(e-e'v, 

B = B{uj,C,,^') is the collision kernel which is positive almost everywhere in its domain 
5^^^ X X R'^, and du) is the normalized measure on the unit sphere S'^~^. 
First of all, we recall the celebrated identity (see [2]) 

(8.2) 4 J mifJ: - fDBdudi'di = |(0 + 0' - 0, - 0;)(/j; - ff)Bdu:di'di 

for any continuous function (j) = Here and below, the integrals are taken over the 

whole domain and we write 0' = = (/>(^*) and 0^ = (piC*)- The identity is a direct 

result of the following symmetry properties of B{uj,C,,C,')'- 

(8.3) B{u;,^,e) = Biu;,e,0 = Biu;,^,,Q. 
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Clearly, the integral in fl8.2p is zero (independent of /) if + 0' = 0^ + 0^. It is well-known 
\2\ that 

(8.4) + 0' = 0, + 0; ifandonlyif 0(0 G span{l, ^i, 6, ■ ■ ■ , ^ 1^1'}- 

Moment closure systems considered here are derived from the Boltzmann equation as 
follows. Let be a positive integer and give n linearly independent continuous functions 
Cfe = Cfc(0 of ,^ G R'^ (fc = 1, 2, ■ ■ ■ ,n). Multiplying (18. ip with Cfc(0 and integrating the 
resulting equations with respect to ,^ G R'^ leads to n equations 

(8.5) dt J Ckfdi + V. ■ y ickfdi = j CkifJi - fnBdudi'di. 

Let Ofc = «fc(a;, t)(A; = 1, 2, ■ ■ ■ , n) be n unknown scalar functions of (x, t). Substituting 

n 

(8.6) / = /(x, t, = exp ( ^ Cfc(0«fc(x, t)) 

k=l 

into (18.51) . we get n first-order PDEs for the n unknown a^. 

Remark 8.1. Traditionally, each is a polynomial of ^ and j c^fd^ is called a moment. 
Here we do not require the c^'s to be polynomials. 

To make clear that the moment closure systems in (18. 5p with (18.61) are of the form (11.11) , 
we write c(0« = J2k=i Ck{C)Oik and introduce the following functions of a G R": 



(8.7) 



exp (c(Oa)c^C, 
Oexp {c{^)a)d^, 



/Qi\ 
Q2 



/Cl\ 

C2 



( exp(c^Q; + c'^a) — exp(c'a + ca))Bdujd^' d^. 



Qia) = 

\qJ \cJ 

Here we have considered 

(8.8) /^/^ = exp(c^Q; + c[a) and //' = exp(cQ; + c'a) 

thanks to the Ansatz in (18.60 . With such an /, we formally have 



drj 

dak 



Thus (18.51) can be rewritten as 

df]a{a 



].9) 



j Ckfd^ and qja^ 
dqja{a) 



at 



dak 



Q{a) 



ijCkfdi. 



Assume that there is a convex open set Q C R" such that the functions in l{8. 7)) are 
well-defined and smooth for a E Q. The existence of such a Q depends on the choice of 
the Cfc's and will not be addressed here. The interested reader is referred to [15] . 

Define 

U = f]a{a) and G = f]aiG)- 
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We show that, for any U E G, there is a unique a E Q such that U = fja{o:). In fact, since 
the Cfc's are hnearly independent, the Hessian matrix 770,0(0;) is symmetric positive definite 
and thereby fi{a) is strictly convex. Then the strictly convex function {fi{a) — a*U) of 
a E Q takes its local minimum at those a satisfying U = f]a{a). Since Q is convex, there is 
at most one such minimum point. Consequently, U = f]a{c() has a global inverse a = a{U) 
for U E G and G is diffeomorphic to the convex open set Q. 
For U eG, set 

ry(f/) := a*{U)U -f]{a{U)). 
We see that the inverse function a{U) is equal to rjij{U). Thus, with 
(8.10) F,{U) := q,a{vu{U)) and Q{U) := Q{vu{U)), 

we arrive at the following system of PDEs: 

dU ^dFJU) , , 

i ■' 

In [15], Levermore showed that ^^(f/) defined above is a strictly convex entropy function 
for (18.111) in the classical sense (I). 

To verify the observed properties (II) and (III), we use (18. 2p and rewrite Q(a) defined 
in (I83D as 



C(a) = 

Notice that 



/ cu + c'l^ - 4 - Ci \ 

+ C'a^ - 4 - C2 



( exp(cv,Q; + c^a) — exp(c'Q; + ca))Bdujd^'d^. 



We set 



and define 



exp(c^a + c'^a) — exp(c'a + ca) 

= / exp[(T(c* + — c' — c)q; + (c' + c)a]da{Ci, + c'^ — c — c)a. 
Jo 

b = b{a, ^, ^*) = / exp[a(c^ + < - c' - c)a + (c' + c)a]rfcr 

^0 



(8.12) aij{a) = (civ^ + c-^ - c ■ - Ci)bB{cj^, + - c^- - Cj)dujd^'d^. 

Thus, C{U) = [aij{rii/{U))]nxn is a symmetric matrix and 

g(t/) = Qivuiu)) = -Ciu)MU). 

Since b and -B are both positive, it follows from f l8.12p that C{U) is non-negative. More- 
over, the null space of C{U) is 

|a G R" : j\ (c(0 + c(e') - c(e.) - c(C))«|'rf^«' = o}, 

which is independent of U. Hence the observation holds for the moment closure systems 
in 1151. 
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9. Discrete Velocity Models 
In this section, we consider discrete velocity models in kinetic theories [7j: 
(9.1) fkt + a{k)-VJk = Qk{U) 

for k = 1,2, ■ ■ ■ ,n. Here fk = fk{x,t) denotes the mass density of gas particles with 
the constant velocity a{k) e R*^ at time t and position x, a{k) ■ = Yl'^j=i(^j{^)'^xj, 
U = (/i, ■ ■ ■ and Qk{U) is the collision term given by 

(9-2) Qk{U) = J2iA'Sf^f)-^Mfkfi), 

ijl 

where the summation is taken over all G {1,2, ■■ ■ ,n} and the coefficients A'lj are 
non-negative constants satisfying 

(9.3) A^h^i = ^i- 

It is not difficult to deduce from these symmetry properties that 

(9.4) Mkiu) = 4'(0fe +<Pl-<P^- <pmfj - fkfi) 

k=l ijkl 

Remark that fl9.3p and (19.41) are analogous to the fundamental properties in (18. 3p and 
(18. 2p of the Boltzmann equation. 

Our aim here is to show that, in the state space 

G:={fk>0: k = l,2,--- ,n} 3U, 

the discrete velocity model (I9.ip - (I9.3I) admits our observation with the strictly convex 
function 

(9.5) r/(f/) = ^/,(log/fc-l). 

fc=i 

The strict convexity of t]{U) is obvious. Since it does not contain any cross-term, t]{U) is 
an entropy function for the diagonal and semilinear system (19. ip . 
To see the properties (II) and (III), we set 

6g = 6g(f/) = / exp[a(log /, + log /, - log fk - log fi) + log fk + log fi]da > 0, 
Jo 

which obviously has the symmetry properties (19. 3p . Then the source terms can be rewrit- 
ten as 

QkiU) = J] Ag6g(log/, + log/,- - log fk - log/0 

(9.6) 

= - 2^ dkm log fm, 
m 

where 

o-kmiU) = — ^ A'^jb'^j — ^ A'll^fi'll^ + ^ Afj^fef™ + 5km ^ A^^hfj 

jl il ij ijl 
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with the Kronecker delta. Therefore, we have 

(9.7) Q{U) = -[akm]nxnVu{U) = -C{U)r]u{U), 

for r]uiU) = (log/i,log/2,--- ,log^)* due to (ED. 

It remains to check the desired properties of C{U) defined in ( 19.71) . Thanks to the 
symmetry properties fl9.3p for both bf- and Af-, it is not difficult to see that akm = cimk, 



that is, C{U) is symmetric. Moreover, let y = (1/1,1/2, 
and (19.61) to obtain 



,yn) G R". We refer to 



ijkl 



Vk - yi) 



> 0. 



Hence, C{U) is non-negative and its null space is 

{1/ G R" : A^,j{yi + yj - y^ - yi) = for all k, l}. 

which is independent of U. Hence, our observation holds for the discrete velocity kinetic 
models constructed in [7j. 

We conclude this paper by writing down the simplified system ( 13. 2p for the discrete ve- 
locity models. To do this, we compute from (19. 5p that rjuiU) = (log /i, log /2, • " " , log /„)*. 
Then the corresponding simplified system reads as 



//A 

/2 







aj{2) 



\ 



//A 

/2 



/log /A 
log/2 



where £ is a constant, symmetric and non-negative definite n x n-matrix. 
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